High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.